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Abstract 

We present preliminary measurements of the CP asymmetries and branching fractions for B — > irir 
and B — > Ktt decays. A total of 347 million BB events collected by the BABAR detector at the 
PEP-II asymmetric-energy e + e~ collider at SLAC are used for these results. We find 
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(5.12 ±0.47 ±0.29) x 10 
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(13.3 ±0.56 ±0.64) x 10 




C^O^O = 


-0.33 ± 0.36 ± 0.08 






-0.019 ±0.088 ±0.014 




A Kn o = 


0.016 ±0.041 ±0.012 



The measured values of 5 T7r and imply that CP conservation in B° — * -k + it decays is excluded 
at the 3.6 a level. From these results we present bounds on the CKM angle a. 
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1 INTRODUCTION 



CP-violating processes are incisive tests of the Cabibbo-Kabayashi-Maskawa (CKM) model of quark 
mixing [1]. Measurements with small theoretical uncertainties in the CKM model provide the 
most effective constraints on physics outside this model. The CKM Unitarity Triangle angle a = 
arg [— ^d^tb/^ud^ubl ^ s measured through the interference of b — > u quark- level decays and B° <-> 
B° mixing. Multiple measurements of a, with different decays, further test the consistency of the 
CKM model. The time-dependent asymmetry in B° — > tt + tt~ is proportional to sin2a in the limit 
that only one amplitude contributes to this decay. However, measurements of an unexpectedly 
large branching fraction for B° — > 7r°7r° compared to that of B ± — > 7r ± 7r° and B° — > tt + tt~ [2] 
indicate that both b — ► u (tree) and b — > d (penguin) amplitudes, with different weak phases, are 
present. Thus, the time-dependent asymmetry is modified to 

\A(At)\ 2 — \A(At)\ 2 
"(At) = ^ At 2 + ^ At 2 = S ™ sin ( Am d A *) ~ Cnn cos (Am d At) 

_ \A\i-\A\ 2 (1) 
nn \A\z + \A\z 
S wn = y/l-C^sm (2a - 2Aa 7r7r ), 

where At is the difference between the proper decay times and Am^ is the f?-meson mixing fre- 
quency. Both the total phase difference Aa n7r and C nn may differ from zero due to a penguin 
contribution to the decay amplitude A. 

The magnitude and relative phase of the penguin contribution to the asymmetry may be un- 
raveled with an isospin relation between the three B — > tttt decays [3]. The amplitudes A tJ for the 
B — ► ■k % tx 3 decays satisfy the relation 

A +0 = -±=A + - + A 00 , (2) 
v2 

with a similar expression for the conjugate amplitudes. The shape of the corresponding isospin 
triangle is determined from measurements of the branching fraction and time-integrated CP asym- 
metry for each B — ► tttt decay. We define the direct CP asymmetry as 

|A 00 | 2 - \A m \ 2 



| A 00|2 + | A 00|2 
_ |i-0|2 _ |^+0|2 W 



|A-°| 2 + |^+°| 2 ' 



From the difference in shape of these triangles for B° and B°, Aa n7r may be determined up to a 
four-fold ambiguity. No CP asymmetry is expected in the AI = 3/2 decay B ± — > ^tt , where no 
penguin amplitudes are present. 

The phenomenology of the B — ► tttt system has been studied in a variety of theoretical frame- 
works and models [4]. Predictions for the relative size and phase of the penguin contribution vary 
considerably, so more precise measurements will help to distinguish among different theoretical 
approaches and add to our understanding of hadronic B decays. 

In addition to the unexpected pattern of decay rates in the B — > tttt system mentioned above, 
the measured rates and direct CP- violating asymmetries in B — > Ktt decays [5, 6, 7] reveal puz- 
zling features that could indicate significant contributions from electroweak (EW) penguins [8, 9]. 
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Various methods have been proposed to isolate the Standard Model contribution to this process in 
order to test for signs of new physics. Sum rules derived from [/-spin symmetry relate the rates and 
asymmetries for the decays B° or B+ to K+tt~, K + ir°, K°ir°, and K°ir + [10], while SU{3) sym- 
metry can be used to make predictions for the Kir system based on hadronic parameters extracted 
from the irir system [8] . 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected in 1999-2006 with the BABAR detector at the PEP-II 
asymmetric-energy B-meson factory at the Stanford Linear Accelerator Center. A total of 347 
million BB pairs were used. The preliminary results presented here supersede the results in three 
prior publications [11]. Roughly 120 million more BB decays have been added, and a number of 
improvements have been introduced to the data analysis, effectively increasing the acceptance for 
the modes containing a neutral pion. 

The BABAR detector is described in detail elsewhere [12]. Charged-particle (track) momenta 
are measured with a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH) inside a 1.5-T superconducting solenoidal magnet. Neutral-cluster (photon) positions 
and energies are measured with an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals. The photon energy resolution is cte/E = {2.3/B( GeV) 1 / 4 © 1.9} %, and the angular res- 
olution from the interaction point is oq = 3.9°/ \/E( GeV). Charged hadrons are identified with a 
detector of internally reflected Cherenkov light (DIRC) and ionization measurements in the tracking 
detectors. The average K-ir separation in the DIRC varies from 12cr at a laboratory momentum 
of 1.5 GeV/c to 2a at 4.5 GeV/c. 

3 ANALYSIS METHOD 

Many elements of the B — > irir measurements are common to the three groups of decay modes 
B° -► h + h'-(h = vr or K), B° -► 7r°7r°, and B ± -► h ± ir°. B candidates (B rcc ) are formed by 
combining two particles, either tracks or ir° candidates. 

3.1 Track and Cluster Selection 

For the B^ — > h^ir and the B° — > h + h'~ samples, we require that each track have an associated 
Cherenkov angle (9c) measured with more than five signal photons detected in the DIRC, where 
the value of 6q must agree with either the pion or kaon particle hypothesis to within 4 a. The 
last requirement efficiently removes events containing high-momentum protons. Electrons are re- 
moved based on energy-loss measurements in the SVT and DCH, and on a comparison of the track 
momentum and the associated energy deposited in the EMC. 

The 7T° candidates are formed from two EMC clusters, one EMC cluster containing two nearby 
photons (merged 7r°), or one EMC cluster and two tracks from a photon conversion to an e + e~ pair 
inside the detector. Previous BABAR results for B° — > 7r°7r° and B^ — > h^ir only included it from 
two EMC clusters; the addition of merged 7r° and converted photons increases the 7r° efficiency by 
10%. Clusters are required to have a transverse energy deposit consistent with a photon, and to 
have an energy E 1 > 0.03 GeV. To reduce the background from random photon combinations, 
the angle # 7 between the photon momentum vector in the 7r° rest frame and the 7r° momentum 
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vector in the laboratory frame is required to satisfy | cos# 7 | < 0.95. The ir° candidates are fitted 
kinematically with their mass constrained to the nominal 7r° mass. 

Photon conversions are selected from pairs of oppositely charged tracks with invariant mass less 
than 30 MeV/c 2 and whose momentum vector points to the beamspot. The conversion point is 
required to lie inside the detector material. Photons from conversions are combined with photons 
from single EMC clusters to form ir° candidates. 

Single EMC clusters containing two photons are selected with the transverse second moment, 
S = Ei x (Ao^) 2 / E, where E{ is energy in each CsI(Tl) crystal and Aaj is the angle between the 
cluster centroid and the crystal. The second moment is used to distinguish merged tt° candidates 
from both single photons and neutral hadrons. 



3.2 Event Selection 

Two kinematic variables are used to separate B decays from the large e + e~ — > qq (q = u,d,s,c) 

background: the beam-energy-substituted mass toes = \J {s/2 + p i • pb) 2 /E 2 — p 2 B , where \/s 

is the total e + e~ center-of-mass (CM) energy, (Ei,pi) is the four-momentum of the initial e + e~ 
system and ps is the -B-candidate momentum, both measured in the laboratory frame, and AE 
= Eb — \/s/2, where Eb is the B-candidate energy in the CM frame. 

Two additional quantities take advantage of the event topology to further separate B decays 
from qq background. The cosine of the angle 9 S between the sphericity axes of the B candidate's 
decay products and that of the remaining tracks and clusters in the event, in the CM frame, 
is peaked at ±1.0 for jet-like qq events, but has a flat distribution for B decays. We require 
|cos# s | < 0.7(0.8) for B° -► ^°7r° {B ± -► /^tt and B° -► h + h'~). For just the B° -► h+h'~ 
sample, we further require that the second Fox-Wolfram moment satisfy i?2 < 0.7 to remove a 
small remaining background from r-pair events. To improve the discrimination against qq events, a 
Fisher discriminant T is formed from the sums ^2 i p% and Y i Pi cos2 @U where pi is the momentum 
and 9i is the angle with respect to the thrust axis of the B candidate, both in the CM frame, of all 
tracks and clusters not used to reconstruct the B meson. 

The number of B decays and the corresponding CP asymmetries are determined in extended un- 
binned maximum likelihood (ML) fits to the variables toes> AE, and J-, plus additional information 
as described below. The likelihood is given by the expression 



M \ N 



exp -J^n, 11 



M 



(4) 



where the product is over the number of events N, the sums are over the event categories M, rii is 
the coefficient for each category as described below, and the probability density function (PDF) V 
describes the distribution of the variables x in terms of parameters a. 



3.3 B° -> 7T+7T- and B° -> K+n~ 

The time-dependent CP asymmetry measurement in B° — > ir + ir~ also uses B-flavor, decay-time, 
and particle-identification information to separate B° — ► it + it~ and B° — > K + ir~ decays and 
measure their CP asymmetries. 

The variables tues and AE are calculated assuming that both tracks are charged pions. B — > -kit 
events are described by a Gaussian distribution for both ttt-es and AE, with resolutions of 2.5 MeV/ c 2 
and 28 MeV, respectively. For each kaon in the final state, the AE peak position is shifted from 
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zero by an amount that depends on the kaon momentum, with an average shift of — 45MeV. We 
require 5.20 < toes < 5.29GeV/c 2 and \AE\ < 0.150 GeV. The large region below the signal in 
toes effectively determines the background shape parameters, while the wide range in AE allows 
us to separate B decays to all four final states in the same fit. 

The Cherenkov angle 9c measured by the DIRC is used to further separate charged kaons and 
pions. The difference between the measured and expected values of 9c, divided by its uncertainty, 
is modeled by a sum of two Gaussian distributions. The 9q PDF is parametrized separately for 
K + , K~ , 7r + , and 7r~ tracks as a function of momentum and polar angle. 

We use a multivariate technique [13] to determine the flavor of the other B meson (-B tag ). 
Separate neural networks are trained to identify primary leptons, kaons, soft pions from D* decays, 
and high-momentum charged particles from B decays. Events are assigned to one of seven mutually 
exclusive tagging categories (including untagged events) based on the estimated average mistag 
probability and the source of the tagging information (Table 1). The quality of tagging is expressed 
in terms of the effective efficiency Q = X^ e fc(l ~ ^ w k) 2 , where e& and Wk are the efficiencies 
and mistag probabilities, respectively, for events tagged in category k. The difference in mistag 
probabilities is given by Aw = Wgo — w^ ■ Table 1 summarizes the tagging performance measured 
in a data sample of fully reconstructed neutral B decays to D^*^(tt + , p + , af) (-Bfl av ). 

The time difference At = Azjfi^c is obtained from the known boost of the e + e~ system 
(/?7 = 0.56) and the measured distance Az along the beam (z) axis between the B iec and B ta , g 
decay vertices. We require |At| < 20 ps and a At < 2.5 ps, where a At is the error on At determined 
separately for each event. The signal At PDF for B° — > tt + it~ is given by 

e -|At|/T f 
/±(At meas ) = ^ {(1 T Aw) 

± (1 - 2w k ) [S n7T sin (Am d At) - cos (Am d At)] } ® R(At mcas - At), 

(5) 

where f£ (/^~) indicates a B° (B°) flavor tag, and where index k is the tagging category. The 
resolution function R(At meSuS — At) for signal candidates is a sum of three Gaussians, identical to 
the one described in Ref. [13], with parameters determined from a fit to the -Bfl av sample (including 
events in all seven tagging categories). The background At distribution is modeled as the sum of 
three Gaussian functions, where the common parameters used to describe the background shape 
for all tagging categories are determined simultaneously with the CP parameters in the maximum 
likelihood fit. 

The ML fit includes eight components: B decays and background with the final states 7r + 7r~, 
K + tt~, K~ir + , and K + K~ . The component coefficients for Kir are parametrized as rix± 7T ^ = 
n Kn (1 T -Ak-k) /2, where Ak-k is the direct CP-violating asymmetry. All other coefficients are the 
product of the fraction of events in each tagging category, taken from -Bfl av events, and the event 
yield. The background PDFs are a threshold function for ttt-es and a polynomial for AE. The 
T PDF is a double Gaussian for the background and an asymmetric Gaussian for the signal. All 
background PDF parameters are allowed to float in the ML fit. 

3.4 B° 7r°7r° 

B° — ► 7T 7r° events are identified with an ML fit to the variables mES; AE, and J 7 . For B° — > 7r°7r°, 
we require itt-es > 5.20 GeV/ c 2 and \AE\ < 0.2 GeV. Tails in the EMC response produce a 
correlation between ttt-es and AE, so a two-dimensional PDF, derived from Monte Carlo (MC) 
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simulation, is used. The PDF for J 7 is a step function, with parameters taken from MC. -Bfl av 
data is used to verify that the MC accurately reproduces the T distribution. The qq background 
PDFs are a threshold function for tties, a polynomial for AE, and a step function for T. All qq 
background PDF parameters are allowed to float in the ML fit. 

The decays B + -► p+vr and B° -► K°ir°(K° -► 7r vr ) add 52 ± 6.8 background events to 
B° — > 7r°7T and are included as an additional fixed component in the ML fit. We model these B 
backgrounds with a two-dimensional PDF to describe toes and AE and with a step function for 
T , all taken from MC simulation. 

The time-integrated CP asymmetry is measured by the P-flavor tagging described previously. 
The fraction of events in each tagging category is also constrained to the fractions determined from 
MC simulation. The PDF coefficient for the B° — > 7r°7r° signal is given by the expression 

»V>7r°,fc = TjfkN^o^l - Sj(l - 2x)(l - 2u/ fc )C 7r o 7r o}, (6) 

where fu is the fraction of events in tagging category k, N n o n o is the number of B decays, x = 
0.184 ± 0.004 [14] is the time-integrated mixing probability, and Sj = +1(— 1) when the B tag is a 
B° 

3.5 B ± -»• ti^ti- and B -> lC±7r 

— > 7r ± 7r° and P 1 * 1 — > K ± 7r° events are identified in an ML fit to the variables tiies, AE, T , 
and 9c- We require B candidates to satisfy m E s > 5.22GeV/c 2 and -0.11 < AE < 0.15 GeV. 
The tighter requirement on AE serves to remove -B-decay backgrounds. The treatment of AE and 
the use of the Cherenkov angle for kaons and pions is identical to that in P° — > ir + -K~. The T 
distribution is also described with a step function, with parameters taken from MC simulation. The 
qq background PDFs are a threshold function for mss, a polynomial for AE, and a step function 
for T . All qq background PDF parameters are floating in the ML fit. 

Several charmless B decays with both a high momentum charged track and a 7r° or photon 
appear as background; significant contributions, as determined with MC, come from the following 
decay modes with the estimated number of events provided in parentheses: B° — > p + ir~ (32 ± 3 
events), B+ -► p+vr (20 ± 2), B -> X sl (5.8 ± 0.6), B ± K ^^ -» 7r°7r°) (5.0 ± 0.3), and 
B° — > p + K~ (3.3 ± 0.5). These events, as well as their measured CP asymmetry, are included as 
an additional fixed component in the ML fit. 

The PDF coefficient for B ± — > Tr 1 ^ and B ± — > K ± 7r° signal is given by the expression 

«i = ^Ni(l - qjAi), (7) 
where Ai is the direct CP asymmetry and qj is the charge of the B candidate. 

4 RESULTS AND SYSTEMATIC UNCERTAINTIES 

Results from the ML fits for the B° — > 7r°7r° and P ± — > /i ± 7r° decay modes are summarized in 
Table 2. Distributions of rriss, AS, and J 7 for B° — > 7r°7r° are in Fig. 1, where a signal-enhanced 
subset of the data is shown. The B ± — > /i 1 * 1 ^ data are shown in Fig. 2. With a large signal in 
both decay modes, we show weighted and background-subtracted plots [15] for the B^ — ► vr^vr 
and 5± i^ ± 7r° signal. The same technique is used to display the qq background as well. 
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The uncertainty in the efficiency for the B° — ► 7r°7r° decay mode is dominated by a 3% systematic 
uncertainty per 7r°, estimated from a study of r — > -ktPv t decays. There is an additional 3.6% 
uncertainty from our knowledge of the EMC resolution function, based on a study of the resolution 
of the 7T° mass and the energy of the photon in e + e~ — > /U + // _ 7 events. Systematic uncertainties 
involving the ML fit are evaluated by varying the PDF parameters and refitting the data. The 
changes in the m,Es and AE signal PDFs are taken from the difference in these quantities in the 
B ± — > h^ir sample between data and MC. The change in the result is taken as the systematic 
error. All sources of systematic error for B — ► 7r°7r° are listed in Table 3. 

The largest uncertainties in the B^ — > h^it decays arise from uncertainty in the hies and AE 
PDFs, our knowledge of the signal T distribution evaluated from a sample of B ± decays, and the 
3% 7T° efficiency uncertainty. The size of a small bias in the ML fit is included as a systematic 
error. The dominant uncertainty on the direct CP asymmetries are taken from the size and error 
in the asymmetry fit in the qq background and the effect of CP violation in the B backgrounds. 
The uncertainties in B^ — > h^ir are summarized in Table 4. 

All results for the B° — > h + h'~ decay modes are listed in Table 5. The correlation coefficient 
between and C n7T is found to be —0.082. The data distributions of m,ES> AE, and T for 
B° -► h + h'- decays are shown in Fig. 3 for signal and Fig. 4 for background with the event- 
weighting technique. The direct CP asymmetry in B° — > K + tt~ is apparent in the distribution 
of AE for B° and B° decays, shown in Fig. 5. We show the distributions of At for signal and 
background decays in Fig. 6. In Fig. 7, we show the distribution of At separately for B° — ► 7r + 7r~ 
events tagged as B° or B°, and the asymmetry a(At). The central values and errors for S nn and 
(7^,,- are shown in Fig. 8, along with confidence-level contours. Our measurement excludes the 
absence of CP violation (S W7T = 0, C W7T = 0) at a confidence level of 0.99970, or 3.6 a. 

Systematic uncertainties for the CP asymmetries Ak-k-, S n7V , and C nw are listed in Table 6. 
For the asymmetry in the K + ir~ mode, we find a background asymmetry of —0.0042 ± 0.0064, 
which is consistent with zero. We therefore take the sum in quadrature of the central value of the 
background asymmetry and its statistical uncertainty as the systematic uncertainty on the signal 
Ak-k to account for possible charge-dependent detector and analysis bias. To further check for 
biases in the fitting technique, we perform a large number of pseudo-experiments where the signal 
events are randomly sampled from simulated MC events, and the background is generated directly 
from the PDFs. We find a bias in Ak-k of 0.002 and include this in the systematic uncertainty. The 
biases on S nn and C n7r in this study are consistent with zero, so we take the sum in quadrature of 
the central value and its uncertainty as the systematic error due to potential bias in the fitter. The 
remaining systematic effects for S n7T and C n7T are dominated by uncertainties in the parameterization 
of 5-flavor tagging and vertexing, and (for C T7T ) in the effect of CP violation on the tag side. 
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5 CONCLUSIONS 



The branching-fraction and CP-asymmetry results described in this paper are: 





-0.53 ±0.14 ±0.02 


C n n = 


-0.16 ±0.11 ±0.03 


Aktt = 


-0.108 ±0.024 ±0.008 


B{B Q -+ vrV) = 


(1.48 ±0.26 ±0.12) x 10 -6 


B(B ± v^vr ) = 


(5.12 ±0.47 ±0.29) x 10~ 6 


B(B ± -► K^ ) = 


(13.3 ±0.56 ±0.64) x 10~ 6 


C^o^o = 


-0.33 ±0.36 ±0.08 




-0.019 ±0.088 ±0.014 


A K7r o = 


0.016 ±0.041 ±0.012. 



Combining these with a branching fraction B(B° — > ir + n~ -) = (5.8±0.4±0.3) x 10~ 6 , also measured 
by BABAR [16], we may evaluate the constraints on both the penguin contribution to a and on the 
CKM angle a itself. Constraints are evaluated by scanning the parameters of interest, |Aa^| = 
| a — a e g | and a, and then calculating the \ 2 f° r the five amplitudes (A +0 , A H , A 00 , A H , A 00 ) 
given our measurements and the isospin-triangle relations [17]. The \ 2 is converted to a confidence 
level (C.L.) as shown in Fig. 9. The upper bound on |Aa^| is 41° at the 90% C.L. Somewhat more 
restrictive new constraints on a are found from measurements of B — > pp and B — > pn decays [18]. 

We have also presented updated preliminary measurements of the branching fraction for K + ir°, 
and the charge asymmetries in K + ir~ and K + ir°. Ignoring color-suppressed tree amplitudes, the 
charge asymmetries in K + 7r~ and K + tt should be equal (see Gronau and Rosner in Ref. [10]), 
which has not been supported by recent data [5]. The values of Ak-k and Aj^^o reported here are 
separated by over two standard deviations. These results could indicate a large color-suppressed 
amplitude, an enhanced electroweak penguin, or possibly new-physics effects [19]. 
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Table 1: Average tagging efficiency e, average mistag fraction w, mistag fraction difference Aw = 
w(B°) — w(B°), and effective tagging efficiency Q for signal events in each tagging category. The 
quantities are measured in the -Bfl av sample. 



Category 






w{%) 




Aw (%) 


Q(%) 


Lepton 


8.67 


± 


0.08 


3.0 


± 


0.3 


-0.2 


± 


0.6 


7.7 


± 


0.1 


Kaon I 


11.0 


± 


0.08 


5.3 


± 


0.4 


-0.6 


± 


0.7 


8.7 


± 


0.2 


Kaon II 


17.2 


± 


0.1 


15.5 


± 


0.4 


-0.4 


± 


0.7 


8.2 


± 


0.2 


KaonPion 


13.8 


± 


0.09 


23.5 


± 


0.5 


-2.4 


± 


0.8 


3.9 


± 


0.1 


Pion 


14.4 


± 


0.09 


33.0 


± 


0.5 


5.2 


± 


0.8 


1.7 


± 


0.1 


Inclusive 


9.6 


± 


0.08 


41.9 


± 


0.6 


4.6 


± 


0.9 


0.25 


± 


0.04 


Untagged 


23.4 


± 


0.12 





















Total Q 30.4 ± 0.3 



Table 2: The results for the B — > ir it and B — > h ir decay modes are summarized. For each 
mode, the number of signal events iVg, total detection efficiency e, branching fraction B, and CP 
asymmetry are given. Errors on B and the asymmetries are statistical and systematic, respectively, 
while errors for iVg are statistical and those for e are purely systematic. 



Mode 


N s 


e{%) 


fi(io- e ) 




Asymmetry 


B^ 


ttV 


140 ± 25 


27.1 ± 1.7 


1.48 ±0.26 ±0.12 


C w o n o 


= -0.33 ±0.36 ±0.08 


B±^ 




572 ± 53 


32.1 ± 1.8 


5.12 ±0.47 ±0.29 


•Ann" ~ 


-0.019 ±0.088 ±0.014 


B±^ 


K** 


1239 ± 52 


26.8 ± 1.3 


13.3 ±0.56 ±0.64 


A Kn o -- 


= 0.016 ±0.041 ±0.012 



Table 3: Systematic uncertainties in the determination of the B — > 7r°7r° branching fraction (left) 
as a percentage change, and the C n o n o asymmetry (right) as an absolute change. 



Source 


A6(7T vr ) 


7r u efficiency 


6.0% 


AE resolution 


3.6% 


ttjes PDF endpoint 


3.1% 


mean of AE and itt-es 


1.9% 


B(B pvr) 


1.8% 


luminosity 


1.1% 


Total 


8.2% 



Source 




tagging 


0.06 


ttzes an d AE 


0.04 


B background asymmetry 


0.03 


Total 


0.08 
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F bin PDFsignai / ( PDF + PDF background ) 

Figure 1: The distributions of (a) mES, (b) AE, and (c) Fisher discriminant T for B — > 
candidates that satisfy an optimized requirement on the signal probability, based on all variables 
except the one being plotted. The projections contain 27%, 30% and 68% of the signal, 20%, 21% 
and 23% of the pir° background, and 2.8%, 0.047% and 4.8% of the continuum background, for 
rriES; AE, and T ', respectively. The PDF projections are shown as a dashed line for qq background, 
a dotted line for B^ — > p^ir and B° — > K°7r°, and a dashed-dotted line for B° — ► 7r°7r° signal. The 
solid line shows the sum of all PDF projections. The PDF projections are scaled by the expected 
fraction of events passing the probability-ratio requirement. Also shown (d) is the ratio of the PDF 
for signal to the PDF for signal plus background comparing data (points) to the components of the 
PDF model. 



18 




F bin F bin 



Figure 2: The distributions and PDF projections of mEs (top), AE (middle), and Fisher discrimi- 
nant T (bottom), for — ► vr^vr (left) and i3 ± — ► K^tt (right) candidates. The main plots show 
the signal data (points) and PDF (line) after event weighting and background subtraction with the 
method described in [15]; the method uses all variables except the one being plotted. The insets 
show the corresponding distributions for background. 
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Table 4: Dominant systematic uncertainties for — > vr^vr and i? 1 * 1 — > X^vr , as percentage 
changes in the branching fractions B (left), and absolute changes in the asymmetries A n ± n o, Ak±k° 
(right). A = (N B0 - N B0 ) I (N B o + N E0 ) . 



Source 


AB(7T ± 7r°) 


Afi(K ± 7T°) 








ttt-es and AS res. 


3.1% 


2.4% 








T PDF 


3.1% 


2.1% 


Source 


A(Anro) 


A(A K7T o) 


7T° efficiency 


3.0% 


3.0% 


tties and AS 


0.007 


0.003 


AS mean 


1.2% 


3.0% 


B background 


0.009 


0.002 


fit bias 


1.1% 


1.8% 


detector asymmetry 


0.008 


0.011 


B background 


0.7% 


0.2% 


Total 


0.014 


0.012 


h^ 1 identification 


0.7% 


0.8% 








Total 


5.6% 


4.8% 









Table 5: The results for the B° — > decay modes are summarized. For each mode, the number 

of signal events iVg and CP asymmetries are shown. Errors are statistical and systematic. 



Mode 


N s 


Asymmetry 


B^ 


7T + 7T 


675 ± 42 


= -0.53 ± 0.14 ± 0.02; C n7T = -0.16 ± 0.11 ± 0.03 


B ^ 




2542 ± 67 


A K7T = -0.108 ± 0.024 ± 0.008 


B ^ 


K+K- 


11 ± 19 





Table 6: Summary of systematic uncertainties on Ak-k-, S nn , and C n7T . The total uncertainty is 
calculated as the sum in quadrature of the individual contributions. 



Source 


Ak-k 


<? 


C 


PDF parameters 


0.0007 


0.0003 


0.0011 


Tagging / Vertexing 




0.0178 


0.0170 


SVT alignment 




0.0100 


0.0022 


Beam spot 




0.0100 


0.0100 


Tag-side interference 




0.0080 


0.0230 


Tgo and Anid 




0.0015 


0.0036 


Potential bias 


0.0078 


0.0051 


0.0036 


Total 


0.0078 


0.0247 


0.0300 
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Figure 3: The background-subtracted distributions of 7Bes (top), AE (middle), and Fisher discrim- 
inant T (bottom) for signal B° — ► it + tt~ (left) and B° — > ET + 7r _ (right) candidates in the data. 
The curves represent the PDFs used in the fit. 
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Figure 4: The signal-subtracted distributions of m-ES (left); AE (middle), and Fisher discriminant 
T (right) for all background h + h'~ candidates in the data. The curves represent the PDFs used in 
the fit. 




-0.1 0.1 



AE (GeV) 

Figure 5: The background-subtracted distribution of AE for signal ET ± 7r =F events, comparing B 1 
(solid) and B° decays (dashed). 
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Figure 6: The background-subtracted distributions of At for signal ir + ir - (left), K tt^ (middle), 
and the signal-subtracted At distribution for background candidates in the data (right). The curves 
represent the PDFs used in the fit. 
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Figure 7: The background-subtracted distributions of At for signal tt + 7t~ events tagged as B° (top) 
or B° (middle), and the asymmetry, defined as A = (N B o — N^ ) / (Nqo + N^ ) (bottom). The 
curves represent the PDFs used in the fit. 
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Figure 8: Central value and errors for S nn and C nn and confidence-level (C.L.) contours. The 
measured value is 3.6 a from the point of no CP violation (S n7T = 0, C nn = 0), converting from the 
C.L. to the units of "two-sided" Gaussian significance. 
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Figure 9: Constraint on the angle Aa„ = a — a e fj (top), expressed as one minus the confidence 
level (C.L.) as a function of |Aa T7r |. We find an upper bound on lAa^l of 41° at the 90% C.L. 
Constraint on the CKM angle a (bottom) expressed as 1 — C.L. The eight peaks correspond to an 
eight-fold ambiguity in the extraction of a; four solutions are from the value and sign of Aa„, 
which is doubled due to the trigonometric reflections between a e g and tt/2 — a e g . Only the isospin- 
triangle relations and the expressions in Eqn. 1 are used in this constraint. The solution at exactly 
a = is excluded at 1 — C.L. = 4.4 x 10~ 5 , not shown in the plot, corresponding to the exclusion of 
*^7r7r = 0; C nn = at 3.6 a. Some of the solutions, and the region around a = 0, can be disfavored 
by other physics information. 
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